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ABSTRACT 

We present an in-depth analysis of stellar activity and its effects on radial velocity (RV) for the M2 
dwarf GJ 176 based on spectra taken over 10 years from the High Resolution Spectrograph on the 
Hobby-Eberly Telescope. These data are supplemented with spectra from previous observations with 
the HIRES and HARPS spectrographs, and V- and R-band photometry taken over 6 years at the Dyer 
and Fairborn observatories. Previous studies of GJ 176 revealed a super-Earth exoplanet in an 8.8-day 
orbit. However, the velocities of this star are also known to be contaminated by activity, particularly 
at the 39-day stellar rotation period. We have examined the magnetic activity of GJ 176 using the 
sodium I D lines, which have been shown to be a sensitive activity tracer in cool stars. In addition 
to rotational modulation, we see evidence of a long-term trend in our Na I D index, which may be 
part of a long-period activity cycle. The sodium index is well correlated with our RVs, and we show 
that this activity trend drives a corresponding slope in RV. Interestingly, the rotation signal remains 
in phase in photometry, but not in the spectral activity indicators. We interpret this phenomenon 
as the result of one or more large spot complexes or active regions which dominate the photometric 
variability, while the spectral indices are driven by the overall magnetic activity across the stellar 
surface. In light of these results, we discuss the potential for correcting activity signals in the RVs of 
M dwarfs. 


1. INTRODUCTION 

M dwarfs are currently highly desirable targets for 
exoplanet surveys, as they allow for detection of 
terrestrial, potentially habitable planets with current 
or upcoming te chnology. In addition to our own 
M dwarf survey (|Endl et al.l I2003L 1200611 , virtually ev¬ 
ery exoplanet search program now dedicates a signifi- 
cant portion of its time allotment surveying M stars (e.g. 

Haghighiponr et ^120101 fBonfils et al.H2ni3HBerta et al.l 

201311 . An outstanding problem surrounding the dis¬ 
covery of the lowest-mass planets is that below RV 
amplitudes of ~ 5 ms“®, stellar activity will cause 
velocity shifts that may imitate or distort the sig¬ 
nal o f an exoplanet d^ieloz et al.l 120011: iHuelamo et all 
20081: iDumnsaue et al.l 1201 2t iRobertson et al.1 12013IJ : 
Santos et al.ll2014ll . This problem is especially acute for 
M dwarfs, since the magnetic activity of old M stars has 
not been studied as thoroughly as for old solar-type stars. 

Nearby M stars are of immense interest for exoplanet 
discovery and characterization. Statistics of Kepler plan¬ 
ets suggest terre strial-size planets should be co mmon 
around M stars ((Dressing fc Gharbonneaul I2013I1 . and 
the observational advantages of M star planets (rela¬ 
tively high RV amplitudes and planet-to-star radius ra¬ 
tios) mean M dwarfs in the solar neighborhood will of¬ 
fer the earliest opportunities to characterize potentially 
Earthlike worlds. As a r esult, upcoming RV instr uments 
such as CARMENES (|Quirrenbach et ^ I2014I1 . HPF 


(iMahadevan et all I2014I1 . and SPIRou ((Artigau et al.l 
l2014f) will focus intensely on these nearby cool stars. 

To date, the most compelling exoplanets found with 
RV around M stars orbit very magnetically quiet 
M dwarfs such as GJ 581 and GJ 667G (although 
even those stars exhibit signif icant RV contributions 
from magnetic activity; see IRobertson et al.l 12014 
IRobertson fc MahadevanlbOUH . On the other hand, new 
dedicated M dwarf RV surveys will target nearby mid- 
late M stars, which tend to be more rapidly rotating and 
magnetically ac tive. In a photom etric activity survey of 
Kepler targets, iBasri et ahl ((20131 1 find that the fraction 
of all M stars more active than the Sun is much higher 
than for hotter stars, exceeding 90% in some temperature 
bins. Some studies have suggested that the RV ampli- 
tudes of stellar signals may be reduced in t he infrared 
((Reiners et al.l[2(Tl0l : iMarchwinski et al.l[2r)T^ . However, 
even low-amplitude noise from activity will be problem¬ 
atic for identifying habitable zone super-Earths. More 
effort must therefore be invested in understanding activ¬ 
ity in late-type stars and its effects on RV. 

GJ 176 is an excellent archetype of a planet-host M 
star with activity levels in between those of very quiet 
stars, such as GJ 581, and the more act ive mid-M 
dwarf s targeted by upcoming experiments. lEndl et al.1 
((200811 initially claimed detection of a 24 Mm planet in a 
10.2-day orbit around the star, but iButler et al.l ((200911 
showed that this orbital solution is inconsistent with the 
Keck/HIRES RV data. The HIRES data did show high 
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Stellar Parameter 

Value 

Reference 

Spectral Type 

M2 

von Braun et al. ('2014') 

V 

9.951 ± 0.012 

Koen et al. (2010) 

K 

5.607 ±0.034 

Koen et al. (2010) 

Parallax 

107.83 ± 2.85 mas 

van Leeuwen (2007) 

Proper Motion 

Pet = 656.85 ± 3.81 mas/yr 



= —1116.20 zb 2.49 mas/yr 

van Leeuwen ('2007) 

Distance 

9.27 ±0.24 pc 


Mass 

0.50 ± 0.02 Mq 

Delfosse et al. (2000) 

Radius 

0.4525 ± 0.0221 Rq 

von Braun et al. (2014) 


3679 ± 77 K 

von Braun et al. (2014) 

Luminosity 

0.0337 ±0.0018 Lq 

von Braun et al. f2014) 

Metallicity ([M/H]) 

0.07 ±0.15 

Schlaufman V Lau^hlin f2010) 


TABLE 1 

Stellar parameters eor GJ 176 


RV variability for the sta r, suggesting an exopl anet might 
still exist in the system. iForveille et al.l (|2009f) presented 
a 2-component solution to their HARPS RV data, in¬ 
cluding a low-mass planet {M s'mi = 8.3 Mq) on an 
8.8-day orbit and an RV signal near 40 days, w hich anal¬ 
yses of photometry (jKiraga fc Steoienl 120071 1. Ha, and 
the Ca II H&K lines showed is the stellar rotation pe¬ 
riod. As a nearby (d = 9.3pc) M dwarf with one known 
super-Earth already, GJ 176 is a very attractive target 
for follow-up searches for additional low-mass planets. 
However, conclusively demonstrating the presence of any 
further planets in this system will require a thorough un¬ 
derstanding and treatment of the stellar magnetic activ- 
ity- 

In this article, we examine the magnetic activity of GJ 
176 using spectral and photometric activity indicators. 
In addition to the previously-observed rotation signal, 
we see evidence of a very long-term activity trend which 
drives a slope in the observed radial velocities. We ex¬ 
plore to what degree this information may be used to 
increase the detection efficiency for planets in the sys¬ 
tem. 


2. STELLAR PROPERTIES 

At a distance of 9.3 parsecs, GJ 176 (M2) is among 
a handful of known exoplanet hosts within 10 parsecs. 
With a V-band magnitude of ^ 10, it is of roughly av¬ 
erage brightness among our M dwarf RV targets. We 
list the complete set of stellar properties for GJ 176 in 
Table [TJ For the mass and metallicity, although more 
recent estimates exist for this star, we use th e photo met¬ 
ric calibration techniques of Delfosse et al.l (1200011 and 
iSchlaufman fc Laughlinl (|2010 ) to remain consistent with 
the stellar characterization presented in our full-sam ple 
activity survey of M dwarfs (|Robertson et al.ll2013all . 

3. DATA 

3.1. Radial Velocity 

We have carried out a dedicated M dwarf radial ve¬ 
locity survey us ing the High Resolution Spectrograph 
IHRS: lTulllll998ll on the Hobby-Eberly Telescope (HET; 
IRamsev et al.l 119^ . Full details on the survey and 
its results, including the obse rving strat e gy and tar- 
get properties can be found in lEndl et ^ (1200311200611 
and iRobertson et'aP (|2013all . The HRS acquires high- 
precision RVs via the iodine cell technique, in which an 


I 2 absorption cell is placed in the light path to superim¬ 
pose thousands of weak, stable I 2 absorption lines over 
the target spectrum. These lines serve as a wavelength 
reference, allowing us to model the stellar Doppler shift 
necessary to produce the observed star-plus-l 2 spectrum. 
Our RVs are extracted u sing the AUSTRAL software 
package (|Endl et al.l 1200011 . 

We have obtained 98 RVs for GJ 176 over 10 years, 
which we list in Table [H We note that earlier veloci- 
ti es based on t he sam e spectra were originally published 
in lEndl et al.l (1200811 . but as we have re-reduced our en¬ 
tire data set with our latest version of AUSTRAL, their 
values may have changed slightly. 

GJ 176 has also been observed intensively by the 
HARPS and HIRES spectrographs. Our RV analysis in 

S uppl emented by HARPS RVs fro m IForveille et al.l 
andIGom es da Si lva et al.l (|2012ll and HIRES RVs 
lutler et al.l (I2009fb 


3.2. Activity-Sensitive Absorption Lines 

Although RV surveys typically monitor the magnetic 
activity of their targets with the Ga II H&K lines, 
the wavelength coverage of HRS in our standard RV 
mode does not extend to those lines. Instead, we 
use the Ha and Na I D lines as activity tracers for 
our M dwarf targets. Ha has been u sed extensively to 
study activity i n low-mass stars fe.g. iKruse et al.l 120101 : 
IBell et ^1201211 . while the use of the Na I D resonance 
feature has more recent ly been shown to be useful as a 
cool-sty actiyity tracer (lAndretta et al.lll99^IDfaz et al.l 
l2007al: iGomes da Silya et al.l 1201 ilk We haye recently 
performed actiyit y analyses of our entire M dwarf sam¬ 
ple using the Ha (|Robertson et al.l[2013aj) and Na I lines, 
and found that while both indices are helpful for studying 
magnetic actiyity, the sodium lines are much more sensi- 
tiveto actiyity that causes RV shifts. In IRobertson et ahl 
dMia), we showed that a magnetic cycle in GJ 328 
causes RV shifts that cause the orbit of its giant planet 
to appear more circular than suggested by the actiyity- 
corrected yelocities. For the case of GJ 176, where it is 
desirable to correct actiyity-related RV signals to iden¬ 
tify low-mass planets, the Na I D lines are particularly 
advantageous. 

We defi ne our sodium indy , /d, according to the def¬ 
inition of iDiaz et al.l (l2007af) . Specifically, we take the 
ratio of the fluxes inside windows centered on each of 
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BJD - 2450000 

RV (ms-i) 

Ij:} {1 A Window) 

Id (0.5 A Window) 

Ina 

^Ca, T 

2935.80776169 

2939.79788804 

2941.98273285 

3254.93830390 

3297.80620165 

-16.44 ±4.66 
-7.46 ±6.05 
-0.46 ±5.38 
-0.92 ±5.67 
-25.22 ±5.62 

0.12488 ±0.00329 
0.12582 ±0.00327 
0.12042 ±0.00315 
0.11368 ±0.00390 
0.12046 ±0.00272 

0.10040 ±0.00236 
0.10184 ±0.00240 
0.09647 ± 0.00215 
0.08739 ± 0.00292 
0.10108 ±0.00224 

0.07098 ± 0.00107 
0.07041 ± 0.00106 
0.07232 ± 0.00115 
0.07385 ± 0.00107 
0.07223 ± 0.00119 

0.02367 ± 0.00074 
0.02368 ± 0.00072 
0.02368 ± 0.00070 
0.02396 ± 0.00080 
0.02392 ± 0.00074 


TABLE 2 

Radial velocities and spectral activity indices for our HET/HRS spectra of GJ 176. The full table will be provided as an 

ONLINE-ONLY SUPPLEMENT TO THE ARTICLE. 


the Na I D lines (Api = 5895.92 A, Ad 2 = 5889.95 A), 
divided by the flux in the nearby pseudocontinuum. We 
measure Id using bo t h 1 A and 0.5 A windows, since 
iGomes da Silva et al.l (1201111 show Id more frequently 
correlates with the .Shk Ca II index when measured with 
a 0.5 A window. In general, we see very little difference 
in the results from between the window sizes, and 
use the 1 A window unless specifically noted otherwise. 
However, there are some minor differences which we will 
discuss in later sections. For reference, we include both 
values of Id alongside their corresponding RVs in Table 
[2 

While we are primarily interested in the sodium fea¬ 
ture, we have also computed /ho, the Ha activity index, 
for each of our HRS spectra. The procedure we use to 
measure Jh^v and its un certainty is documented fully in 
[Robertson et al.l (|2013all ; like Id , Iro is simply the ratio 
of the flux in a 1.6 A window centered on the Ha line to 
the nearby pseudocontinuum. 

As a control quantity, we measure the flux index 
/cai for the Ga I line at A = 6572.795 A. This line is 
not sensitive to stellar activity, and should therefore re¬ 
main roughly constant. All of our absorption-line indices 
may be found alongside the RVs in Table [2] 

To enable comparison, we have extracted /d and 
Ilia values (where applicab le) from the HARPS and 
HIRES spec tr a disc ussed in iForveille et al.l (I2009D and 
iBiitler et al.l (|2009[ 1 via the escQ and Keck pub¬ 
lic archives, respectively. We also consider activ- 
ity indices from the HARP S M dwarf activity survey 
(|Gomes da Silva et al.l[2Ml) . 

Unfortunately, the Na I D feature lies on the edge of a 
spectral order on HIRES, so we cannot obtain Id for the 
HIRES data. At the edge of the order, the blaze function 
is at a minimum, so any recovered Id values would be 
very low S/N, and the pseudocontinuum against which 
we measure the index would be even more unreliable. 
Furthermore, barycentric velocity shifts frequently cause 
one or both Na I lines to fall off the GGD completely, 
making Id measurements impossible. 

3.3. Photometry 

Prompted by the presence of significant stellar activity 
in its RVs, we have monitored GJ 176 for photometric 
variability with the Tennessee State University (TSU) 
automated Celestron C-14 telescope. The telescope was 
equipped with an SBIG STL-IOOIE GGD camera observ¬ 
ing through Cousins R and Johnson V filters. During the 

^ Based on data obtained from the ESO Science Archive Facility 
under request number 103236. 


BJD - 2450000 

dR 

5854.8932 

-2.43229 ± .00044 

5856.8425 

-2.43041 ± .00061 

5857.0246 

-2.42685 ± .00087 

5858.0320 

-2.42596 ± .00026 

5859.9542 

-2.42502 ± .00063 


TABLE 3 

Differential H-band photometh.y of GJ 176. Magnitudes 

ARE reported RELATIVE TO THE MEAN OF FI VE CONSTANT STARS 
IN THE SAME GGD FIELD (SEE SECTION |3.3ll . ThE FULL TABLE 
WILL BE PROVIDED AS AN ONLINE-ONLY SUPPLEMENT TO THE 
ARTICLE. 


BJD - 2450000 

dy 

4355.7962 

-2.02130 ±.00110 

4355.8337 

-2.02330 ± .00090 

4355.8823 

-2.01750 ± .00100 

4355.9249 

-2.02290 ± .00140 

4356.9144 

-2.01950 ± .00350 


TABLE 4 

Differential V-band photometry of GJ 176. Magnitudes 

ARE reported RELATIVE TO THE MEAN OF FI VE CONSTANT STARS 
IN THE SAME GGD FIELD (SEE SECTION |3.3ll . ThE FULL TABLE 
WILL BE PROVIDED AS AN ONLINE-ONLY SUPPLEMENT TO THE 
ARTICLE. 


2007-2008 observing season, the C14 was located at Van¬ 
derbilt Univerity’s Dyer Observatory in Nashville, Ten¬ 
nessee, where it collected 426 observations of GJ 176 in 
the Johnson V pass band. In 2010, the telescope was 
relocated to TSU’s automated telescope observing site 
at Fairborn Obser vatory in the Patagonia Mou ntains of 
southern Arizona (lEaton. Henry, fc Fekell[2003D . There, 
it collected 42 Johnson V observations in the 2010-2011 
observing season, 137 V and 87 Cousins R observations 
in 2011-2012, and 153 R observations during 2012-2013. 

The i?-band and U-band differential magnitudes are 
listed in Tables [3] and SI respectively, and plotted as a 
function of Julian date in Figure [TJ Fach observation 
has been corrected for bias, flat-field, differential extinc¬ 
tion, and pier-side offset. The differential magnitudes 
are computed against the mean of 4 constant compari¬ 
son stars identified in the same CCD field of view. Fach 
differential magnitude is the mean of 4 to 10 successive 
frames taken on a given night. Period analyses of the 
i?-band and U-band observations are shown in Figures!^ 
and[3l 
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Photometric Phase 

(b) 


Fig. 1.— a. Differential V- and i?-band photometry of GJ 176. The top three panels show the data taken in the V band from Dyer (top panel) and Fairborn (panels 2 3) 

Observatories. The bottom two panels show the R band data from Fairborn Observatory. The middle of the each observing season (opposition) occurs roughly around the center of 
the X axis. The season mean brightness is indicated by the dotted line in each panel, b. Data from a, normalized and combined into a single data set [top). The periodogram for the 
complete set is given in the middle panel. In the bottom panel, we fold the data to the rotation period, with the best-fit sinusoidal model to the data shown as a solid curve. 
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2455900 2456000 2456100 2456200 2456300 2456400 

Barycentric Julian Dale 


(b) 

Fig. 2.— a. Generalized Lomb-Scargle periodograms for our R- 
band photometry. The top panel gives the periodogram of the 
original data, while the middle three panels give the power spec¬ 
trum after subtracting a sinusoidal fit to the 40-day stellar rotation 
period (residual 1), rotation plus the 112-day periodicity (residual 
2), and rotation, the 108-day signal, and the 47-day signal (resid¬ 
ual 3). The fifth panel shows the window function, or the peri¬ 
odogram of our time sampling. The dash-dotted lines indicate the 
power required for a false alarm probability of 0.01 according to 
our bootstrap FAP estimate, b. Phase plots of the periodic signals 
identified in (a), and the residuals around the 3-signal fit (bottom). 
Our sinusoidal fits to the data are given as black lines. 


4. STELLAR ACTIVITY 

We begin our stellar activity analysis by searching for 
periodic behavior in the photometry and/or the spectral 
line indices. In order to examine each data set in its en¬ 
tirety, we use t he generalized Lomb-Scargle pe riodogram 
as described in [Zechmeister fc Kursted (120091 1. Adapted 
from the Fourier transform power spectrum, the Lomb- 
Scargle periodogram excels at identifying periodic behav¬ 
ior in data taken with uneven time sampling. The gener¬ 
alized version from Zechmeister & Kiirster additionally 
allows for individually-weighted data points and floating 
means. 

Because the power of a peak in a Lomb-Scargle peri- 


Period (days) 

FAP 

FAP 


Photometric Signals 



R band 

V band 

39.5 


< 10-4 

< 10-4 

112 


< 10-4 


47.4 


< 10-4 


18 



< 10-4 

70 



< 10-4 



7d Signals 

1 A window 

0.5 A window 

73 


0.0019 

< 10-4 

112 


0.0055 



TABLE 5 

False alarm probabilities (FAPs) for periodic signals 

OBSERVED IN OUR DATA. FAPS ARE B ASED ON 10 * BOOTSTRAP 
RESAMPLING TRIALS, AS DESCRIBED IN IKURSTER ET AlI II1997I1 . 


odogram is related to, but not directly indicative of its 
statistical significance, we estimate false-alarm probabil¬ 
ities (FAPs) for our candidate sign al detections using th e 
bootstrap resampling technique of lKfirster et al.l (1199711 . 
The method retains the time stamps of the original data 
set, while drawing at random (with replacement) a value 
for each time from the set of observed values. The peri¬ 
odogram is computed for a large set of such “fake” data 
sets, and the FAP is taken as the number of resampled 
periodograms with a peak at any period stronger than 
that of the candidate signal. For each signal discussed be¬ 
low, we have computed a FAP in this manner and listed 
it in Table [5l 

We show the results of our R- and P-band photometry 
in Figures [2] and [SI respectively. In both bandpasses, we 
see a strong periodogram peak near 39.5 days, which has 
previ ously been shown to be the rotation period o f GJ 
176 (jKiraga fc Steoieni l2007t iForveille et al.l 1200^ . In 
each bandpass, we have fit a sinusoid of the form F{t) = 
Fq + Asin(wt + (j)), where Fq is the mean value, A the 


amplitude, uj = ^ the angular frequency, and (j) is the 
phase. In i?, the rotation signal has a period of 39.61 ± 
0.07 days and an amplitude of 11.7 ± 0.5 millimags. In 
V, the period decreases to 39.44 ± 0.01 days, while the 
amplitude is 12.0 ± 0.4 millimags. 

In order to obtain the most precise estimate of the 
rotation period possible, we normalized the R- and V- 
band photometry to a mean value of 1.0 and com bined 


them into a single data set, which we show in Figure 1(b) 
The combined photometry yields a single coherent signal, 
for which we derive a period of 39.457 ± 0.011 days. A 
phase plot of the combined photometry at th e rota tion 


period is shown in the bottom panel of Figure 1(b) 
Upon fitting and subtracting the 39-day signal, we see 
a number of additional signals in the residual photome¬ 
try. The i?-band residuals show peaks near 112 days and 
47 days, while the U-band data show an 18-day harmonic 
of the rotation period. After removing the 18-day signal, 
the U-band data show a third periodogram p eak at 70 
days . We include fits to these periods in Figures 2(b) and 


3(b) When including the 47-day signal in our model of 
the A-band data, the period of the longest signal shifts 
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Fig. 3.— a. Generalized Lomb-Scargle periodograms for our V- 
band photometry. The top panel gives the periodogram of the 
original data, while the middle three panels give the power spec¬ 
trum after subtracting a sinusoidal fit to the 40-day stellar rotation 
period (residual 1), rotation plus its 18-day harmonic (residual 2), 
and rotation, the 18-day harmonic, and the 70-day signal (resid¬ 
ual 3). The fifth panel shows the window function, or the peri¬ 
odogram of our time sampling. The dash-dotted lines indicate the 
power required for a false alarm probability of 0.01 according to 
our bootstrap FAP estimate, b. Phase plots of the periodic signals 
identified in (a), and the residuals around the 3-signal fit (bottom). 
Our sinusoidal fits to the data are given as black lines. 


slightly, to 108 days. However, since the periodicity ap¬ 
pears consistently at 112 days in multiple indicators, we 
will continue to refer to it as a 112-day signal. We note 
that our residuals to a 3-sinusoid fit to the V-band pho¬ 
tometry contain some marginally significant power be¬ 
tween 10 and 13 days. Attempting to include a fourth 
sinusoid in our model does not significantly improve the 
fit, and there is no clearly preferred period in this range. 
We therefore conclude that this excess power is likely 
noise. 

Of the four periods identified in the residual photome¬ 
try, two can be easily explained. We interpret the 18-day 
harmonic as evidence that spots or spot complexes oc¬ 


casionally appear at opposing longitudes on the stellar 
surface, creating periodicity at half the rotation period. 
The 47-day signal is close to both the rotation period 
and its 1-year alias (44.2 days). It is possible that differ¬ 
ential rotation creates periodicities near the fundamen¬ 
tal period, as suggest ed for the false-pos itive exoplanet 
signal for HD 41248 (iSantos et al.ll20l^ . The 70 and 
112-day signals, on the other hand, are somewhat sur¬ 
prising. Rotating starspots are expected to create signals 
at the rot ation period and it s integer ratios (i.e. Prot/2, 
Prot/^- ■ ■, iBoisse et al.ll201lll . so signals longer than the 
rotation period but shorter than is typical for magnetic 
cycles are puzzling. We will discuss these periodicities 
further in later sections. 

The periodogram for GJ 176’s Id series (Figure 4(a) I 
shows a number of significant peaks. We discard peaks 
near 1 and 2 years as aliases caused by our observing 
cadence; fits to those periods reveal large phase gaps, as 
is typical of sampling-related aliases. More compelling is 
the presence of significant power at longer periods. While 
it is possible to fit a complete cycle with P = 1200 day 
sinusoid, the RMS around the fit is essentially equal to 
that of a straight-line fit, causing us to adopt a linear 
trend with slope (5± 1) x 10“® day“^ as our best model 
to the data. 

Regardless of the model adopted for the long-term be¬ 
havior of the Na I feature, an additional periodogram 
peak remains in the residuals at P = 73 days. A boot¬ 
strap FAP test produced no false positives in 10^ iter¬ 
ations, leading us to conclude the signal is statistically 
significant. Furthermore, the period of the peak is an ex¬ 
cellent match to the 70-day signal found in our V-band 
photometry. Since 73 days is close to twice the 39-day 
rotation period, we speculate this signal may be related 
to the stellar rotation, although again, such behavior is 
not predicted by simple spot models. 

In analyzing the residual signals in Id , we find a small 
dependence on the size of the windows used to com¬ 
pute the index. While the 73-day peak appears in the 
periodogram regardless of the window used, the S/N 


is slightly higher when using 0.5 A windows. On the 
other hand, the 1 A windows show some power near the 
112-day period observed in the i?-band photometry, al¬ 
though this detection is less significant. Because of this 
slight discrepancy, we sh ow th e Id periodograms from 
both windows in Figure |4(a)| Because of the higher 
S/N, though, we compute our fit to the 73-day sig¬ 
nal using the 0.5 A windows. When including a sine 
curve model alongside the linear trend, we find a fit with 
P = 73.4 ± 0.4d and A = 0.0052 ± 0.0001 (4%). 

Unlike the photometry and the sodium index, 
/net appears to be devoid of coherent periodic behav¬ 
ior. We show the /ho time series an d its periodogram 
in Figure El It is well established le.g. iCincunegui et al.l 
120071 : iGomes da Silva et al.ll2011h that Ha often does not 


correlate with other spectral activity tracers, so the ab¬ 
sence of the long-term slope or the rotation signal does 
not imply a fault in o ur data or analysis. However, since 
iForveille et al.l (1200911 recovered the 39-day rotation pe¬ 
riod in Ha, it is important that we point out our overall 
data set does not show the same signal. 

It is curious that Id should show a multiple of the 
rotation period instead of the period itself. It is also po- 
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Fig. 4. — a. Generalized Lomb-Scargle periodograms for /q, our Na I D index. The top panel shows the periodogram of the original 
time series. While the power spectrum shown is for the computation of Id using 1 A line windows, we note that the periodogram of the 
0.5 A index shows no significant difference. The two subsequent panels show the residual periodogram after fitting and removing a linear 
trend from the original data. Here, we include the periodograms from each of the two Id indices to show that the 73-day period is the 
only peak that remains constant in each. The residual periodogram to a trend-plus-sinusoid fit of the 0.5 A Id and the window function 
are included in the bottom panels. The horizontal lines indicate the power required for a false alarm probability of 0.5 (dot), 0.1 (dash), 
and 0.01 (dash-dot) according to our bootstrap FAP estimate, b. Top: Time-series Id values from our RV survey. Our model of the linear 
trend plus a 73-day sinusoid is given as a solid red line. Below each Id value is the residual around the 2-signal fit and the corresponding 
Ca I index, which is an activity-insensitive line used as a control. Bottom: Residual Id values after removing the linear slope, folded to the 
73-day period identified in the periodograms. We note that while the 73-day peak appears in both Id indices, we show the 0.5 A values 
here 
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Fig. 5. — a. Generalized Lomb-Scargle periodograms for Ina^ our Ha index. The dash-dotted line indicates the power required for a 
false alarm probability of 0.01 according to our bootstrap FAP estimate, b. Our time-series /ho data. We see no periodic signals over the 
entire set of our observations in Ha. 


tentially disconcerting that the rotation period does not 
appea r to be present in /hq, either, since iForveille et al.l 
(j2009f l detected it in both Ca II H&K and Ha. Exam¬ 
ining lua and Id from the Forveille spectra, we confirm 
that both Ha and Na I D show a periodic signal at the 
rotation period. Furthermore, we note that the 112-day 
periodicity we observe in R appears in the HAR PS spec¬ 
tra in Ha and Ca II H&K (jForveille et al.l[200^ Fig. 6). 

The HARPS data cover a relatively short time span 
compared to our own, and it is possible that the rotation 
signal does not stay constant in phase for the absorption¬ 


line indices. We have attempted to verify this hypoth¬ 
esis by restricting the analysis of our data to a subset 
of dense time sampling from September 2008 to March 
2010. While this time frame is still considerably longer 
than a single rotation, it is the shortest time over which 
we have enough data to perform meaningful frequency 
analysis. In order to properly evaluate short-period be¬ 
havior, we have also removed the long-term trend from 
the Id data. 

In Figureini we show periodograms for I^a and Id over 
this abbreviated time period. Here again, we use the 0.5 

































































































Fig. 6.— We have examined the subset of our spectra spanning September 2008 and March 2010 to search for evidence of the stellar 
rotation signal in /hq and Jd. Periodograms of both /hck (a) and (6) show clear peaks near the 39-day rotation period. The dash-dotted 
lines indicate the power required for a false alarm probability of 0.01 according to our bootstrap FAP estimate. We include phase plots of 
the rotation signal in c and d. The presence of the rotation signal in this truncated data (but not in the complete data) suggests the signal 
does not remain in phase over the duration of our observations. 


A window for the Na I index because of its slightly higher 
S/N. Interestingly, both indices show distinct peaks near 
the 39-day rotation period, each at a p ower level sim¬ 
ilar to the absorption-line detections of iForveille et alJ 
(j2009l l. We fit sinusoids to these peaks, which we in¬ 
clude in Figure 0 finding best-fit solutions at P = 34.8 
days (Ha) and P = 37.5 days (Na I). We note, however, 
that there are a number of periods between 35 and 40 
days that are consistent with these data, since the fit 
is not well constrained with only 40 points. Based on 
these detections, we conclude that the rotation period is 
present in our absorption-line indices, but does not re¬ 
main constant in phase, leading to non-detections when 
examining the data in its entirety. 

While the baseline of the HARPS observations make 
assessing the presence of the long-term I-q trend diffi¬ 
cult to assess conclusively, we note that the HARPS 
Id values are consistent with the slope we observe with 
HR S. In Figure [3 we show our measurements of Id from 
the IForveille et al.l (1200911 spectra alo ngside those mea- 
sured bv IGomes da Silva et al.l (1201211 . We have scaled 
our values by a factor of 1.738 to account for a difference 
in normalization between our measurements and those 
of the HARPS team, and exclude two spectra which ap¬ 


pear to show flare eventfl While the observed slope 
((8 ±4) X 10“® day“^) is only significant at the 2a level, 
it is entirely consistent with the measured HRS trend. 

5. RADIAL VELOCITY 

GJ 176 has been monitored extensively by HRS, 
HARPS, and HIRES, resulting in a robust RV data set 
for the star. In Figure [51 we show all 196 available ve¬ 
locities for the star. 

5.1. The HRS and HIRES data 

Because the HRS and HIRES RV series have similar 
baselines, precisions, and results, we will discuss them 
together before comparing/contrasting with the HARPS 
RVs. 

As GJ 176 is relatively nearby {d = 9.3 pc), its secular 
acceleration of ~ 0.4 m s“^ yr“^ has become significant 
over our 10-year observational baseline. Upon subtract- 
ing this acceleration (following iZechmeister et all l2009fl 
from our RV^, we see that a positive linear trend still 

2 BJD = 2453367.7, 24 53814.5 

® IForveille et al.l II2009I') subtracted secu lar acceleration from the 
HARPS velocities, but lSutler et all II2009I) did not subtract it from 
the HIRES data. We have therefore removed the secular accelera¬ 
tion from the HIRES velocities. 
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Fig. 7.— /d from the HARPS spectra of IForveille et aP Il200fll . our measurements, black) and IGomes da Silva et ahl II2012I . blue). The 
measured linear trend (red) is consistent with that seen in the HRS data (cyan). The zero-point offset between the HARPS and HRS slopes 
is caused by uncertainty in the normalization of Iy> used by HARPS, and has no physical significance. 


remains in our v elocities. Plotting RV as a function of 
/d (Figure 9(a) I, we see from the correlation between 
the two variables that this slope appears to be the RV 
signature of the stellar activity trend. Removing the lin¬ 
ear fit from the relation eliminates the trend from the 
RVs, suggesting the acceleration is most likely caused by 
a long-term magnetic cycle rather than an unseen binary 
c ompanion. 

iButler et al.l (|2009D indicated the possibility of a linear 
trend in their Keck/HIRES velocities of GJ 176 which, 
if confirmed, would add confidence to our own detection. 
Examining all the available velocities (Figure [5]), while 
the data show considerable scatter-as expected from the 
planetary and rotation signals-a positive linear trend is 
indeed present over the entire RV series. Specifically, we 
find a linear least squares fit of 


RV(m s"i) = -1.1(0.8) -b 0.0013(0.0005)1 (1) 

Here, t = BJD — 2 456 300. This fit allows for a zero- 
point offset between the two data sets. We obtain a 
Pearson correlation coefficient r = 0.2 for the RV trend 
which, for a set of 141 total RVs, gives a probability 
P{r) = 0.009 of a flat slope. Fitting to the HRS or 
HIRES sets individually yields fits consistent to within 
the uncertainties with Equation [TJ We include all three 
fits in Figure [51 

Although the high statistical significance of the 
/d trend and its correlation with our HET/HRS veloc¬ 
ities convinces us the overall RV set should contain a 
slope, because of the large scatter in the RVs, the shallow 
slope derived from the combined data is only a ^ 2.5cr 
detection. In addition to causing large deviations from 
the fit, the scatter increases the uncertainties in the zero- 
point offsets between the RV sets, which can have a large 


effect on the linear fit. In an attempt to remedy this 
problem, we attempted computing Equation |T] using the 
residual RVs around a fit to the planet. In addition to 
reducing the velocity scatter, fitting and removing the 
exoplanet signal offers a more reliable way to determine 
the zero-point offsets, since the presence of a slope in the 
RVs could otherwise lead to incorrect results. Reassur¬ 
ingly, our fit to the slope is identical regardless of whether 
or not we subtract the planet. This result, coupled with 
the fact that removing the slope from the RVs increases 
the significance of previously-known signals (see below), 
lends additional evidence to the veracity of Equation [T] 

Because we do not have Id values from HIRES, we 
cannot verify the RV-/d correlation for the HIRES data. 
However, given the agreement of the HIRES slope with 
that seen in the HRS data (which we know is correlated 
with activity), we strongly suspect the observed trend is 
due to stellar magnetic activity. With activity-induced 
RV contributions from stellar rotation and a long-term 
magnetic cycle, GJ 176 joins a Gen B among stars with 
multiple stellar RV signals and (candidate) planets. 

While the rotation and exoplanet signals add scat¬ 
ter and inflate the uncertainties on the RV-/d relation 
(or its manifestation as an RV trend), it is neverthe¬ 
less tempting to attempt an activity correction for this 
star to ev aluate how the remaining RV signals respond. 
In Figure |9(b)| we show periodograms of the combined 
HRS-bHIRES RVs before and after subtracting Equation 
dl The power of the peak corresponding to the planet 
{P = 8.78 days) increases from 27 to 28 after subtract¬ 
ing the trend which, assuming the false alarm probability 
for power Z scales as \fZ^ translates to a reduction 
in FAP by a factor of ^ 3. 
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Fig. 8.— All available RVs for GJ 176. The solid line shows our best linear fit to the combined HRS+HIRES data set, while the dashed 
pink and dot-dashed green lines show fits to the HRS and HIRES RVs, respectively. Secular acceleration has been subtracted from the 
velocities, and we have fit and subtracted a constant offset for each data set. 


5.2. The HARPS data 

At first glance, our results for both /d and RV in GJ 
176 appear to contradict those from the HARPS sur¬ 
vey. However, looking at all the available data for the 
star with careful considerations for instrumental preci¬ 
sion and time baselines shows the findings presented here 
a re not necessarily incompatib le with previous studies. 

iGomes da Silva et al.l (1201211 list GJ 176 as inactive in 
Id because it does not pass their variability test, but 
this result can be attributed mostly to their sampling 
and analysis methods. Their data for the star consists of 
just 7 nights of observations, each binned into a nightly 
average. The approach is designed to search for long¬ 
term variability and RV correlation. We therefore do 
not expect the HARPS data to show rotationally-induced 
modulation, as such behavior is intentionally ignored. 

Interestingly, although the combined RV set shows the 
39-day stellar rotation signal at statistically significant 
power, when examining each individual data set we find 
that the HARPS velocities are the only ones which con¬ 
tain any periodogram power at the rotation period. We 
attribute this feature to the fact that the HARPS obser¬ 
vations have a much more frequent cadence, and there¬ 
fore have a higher probability of resolving a rotation sig¬ 
nal before it changes phase. Additionally, since HARPS 
has a much higher resolving power (i? ^ 110 000) than 
HRS and HIRES {R - 50 000-60000), it may be more 
sensitive to spot-induced modulation of the stellar line 
profiles, leading to a more robust detection of the rota¬ 
tion period. The relatively short time baseline prevents 
phase shifts from degrading the rotation signal, and the 
relatively high precision causes the signal to persist in 
the combined RV set since generalized Lomb-Scargle pe- 
riodograms and minimization algorithms give added 


weight to the HARPS RVs. 

It is important to note that the absence of the rotation 
signal from the HRS/HIRES RVs cannot be attributed 
to RV precision alone. Performing a 2-signal RV fit to the 
HARPS velocities alone, we find an RV amplitude of 4.4 
ms“^ for the rotation signal, compared to 4.0 ms“^ for 
the planet signal. As evidenced by the strong detection 
of the planet in the combined HRS/HIRES data, if the 
rotation signal remained constant in phase and ampli¬ 
tude, it should have been easily recovered. 

We are unable to confidently determine the presence of 
the long-term RV slop e for the HARPS RVs. In b o th the 
iForveille et al.l (120091 1 and IGomes da Silva et al.l (|2012[ 1 
RV sets, the raw velocities actually show a negative slope, 
seemingly contradicting the HRS/HIRES slopes. How¬ 
ever, the residuals to a planet-plus-rotation model ap¬ 
pear flat. Evidently, the planet and (especially) the ro¬ 
tation signal dominate the RVs during this period, and 
have a strong effect on any long-term trends, or lack 
thereof. Thus, although the HARPS Id values are con¬ 
sistent with the positive activity trend, we are unable to 
confirm whether this consistency extends to the RVs due 
to the mitigating factors of the short time baseline and 
the periodic signals. 

The HARPS data for GJ 176 create a conundrum for 
this analysis. The HARPS RVs are the most precise of 
the three data sets, and the dense time sampling is ideal 
for resolving the signals of the planet and the stellar rota¬ 
tion. On the other hand, the relatively short time base¬ 
line and the unusually large RV contribution from stellar 
rotation make it difficult to compute a unified model of 
the planet, the stellar rotation, and the long-term trend 
for the combined HRSj-HIRES-l-HARPS RVs. Even if 
the rotation-induced RVs were present throughout the 
observational baseline, phase shifts of the rotation sig- 
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nal would likely prevent a “global” 2- or 3-signal model. 
Since we are primarily interested in the magnetic activity 
of GJ 176 in this study, we have elected not to pursue a 
combined RV fit for this reason. 

6. DISCUSSION 

Counting the stellar rotation period, its first harmonic, 
the long-term trend and the two intermediate-period sig¬ 
nals, GJ 176 exhibits at least five distinct stellar activity 
signals. Two (rotation and the trend) have already been 
observed to manifest in RV, and it is reasonable to expect 
that at the RV amplitudes (< 1 ms“^) of potential addi¬ 
tional terrestrial planets in the system, the other activity- 
related periodicities will appear. Thus, while GJ 176 is 
certainly an interesting candidate system for discover¬ 
ing and characterizing low-mass exoplanets, its activity 
makes the RV interpretation especially difficult. 

With multiple activity-induced RV signatures, GJ 176 
has a reputation as a highly active M dwarf. However, 
it is important to note that with a mean [Tna/Aboi] of 
-3.84, GJ 176 has an overall stellar activity level that is 
es sentially average for its mass (according to the relation 
in IRobertson et al1l2013all . Its rotation period of 39 days 
is likewise ordinary for an old M star, as opposed to the 
shorter rotation periods of more active M dwarfs. Rather, 
what is remarkable about GJ 176 is the fact that its ma¬ 
jor spots have survived for at least the 6 years covered 
by our photometric observations, resulting in the coher¬ 
ent photometric signal. Our results seem to suggest that 
atypically strong magnetic fields (such as would produce 
abnormally large mean activity levels) are not required 
to preserve such spots. 

It is possible that individual spots need not survive 
for extended periods in order to create the coherent pho¬ 
tometric signal we observe. On the Sun, spots prefer¬ 
entially appear at so-called “active longitudes,” where 
increased magnetic activity in a localized region causes 
spots to manifest re peatedly (e.g. lBerdvugina fc UsoskinI 
l2003t llvanovi [200711 . Active longitudes rotate in phase 
with the stellar rotation (modulo differential rotation), 
and coul d explain a pe rsistent coherent starspot signal. 
However. Ilvanovi (120071 1 finds that prominent solar active 
longitudes tend to survive for about 20 rotations or less. 
If the ~6-year photometric signal of GJ 176 is caused 
by a single active longitude, it has survived at least 50 
rotations, suggesting a significant departure from solar 
behavior. More detailed work is required to understand 
how such long-lasting active regions might be maintained 
on low-mass stars. 

The physical origins of the ^73- and ~112-day activity 
signals are of some interest, as they are not simple har¬ 
monics of the rotation period. However, they are close 
to 2 and 3 times the rotation, so it is possible the signals 
are physically related to rotation. 

A study of starspo ts on the M4 dwarf GJ 1243 
([Davenport et al.l 1201411 offers insight towards a poten¬ 
tial explanation. The Kepler lightcurve of GJ 1243 
shows sinusoidal variations induced by stellar rotation 
remaining in phase over more than 4 years, while a sec¬ 
ondary “shoulder” in the lightcurve changes in phase 
on ^ 100-day timescales. These features may be ex¬ 
plained by one or more major spots (or perhaps an active 
region/longitude) persisting on the stellar surface over 
many rotations, while small, short-lived spots change 


over the 100-day timescale. 

In the case of GJ 176, we propose a similar scenario. 
One or more large, highly-spotted active regions must 
survive over many years, dominating the photometric 
variability and leading to the observed in-phase rotation 
signal in the photometry. In addition to the active re¬ 
gions, many “minor” spots can appear at variable times 
and latitudes outside the active regions, also tracing the 
stellar rotation, but not in phase over long timescales. 
The 73- and 112-day (or 2-3 rotation) periods would 
then represent typical minor spot lifetimes. 

So why does the stellar rotation signal remain in phase 
in RV and photometry, but not in the absorption-line 
fluxes? We suggest this is due to the fact that the 
photometry traces starspots (and thus rotation) directly, 
whereas the absorption-line indices trace emission due 
to the magnetic activity producing the spots. If the 
magnetic field across the stellar surface is changing on 
the timescale of the individual spot lifetimes, then those 
changes may dominate over localized active longitudes in 
the spectral activity indices. 

Our activity-RV analysis confirms the planetary nature 
of the 8.8-day periodicity. The planet’s period is not an 
integer-ratio harmonic of the rotation period, remains in 
phase, and does not appear in any of our activity tracers, 
leaving no reasonable suspicion that the signal is pro¬ 
duced by activity. On the other hand, the phase-shifting 
of the 39-day signal in RV, /hq and Id ensures it is in 
fact the rotation period, and not produced by magnetic 
interaction between GJ 176 and a second, more distan t 
planet (a possibility mentioned bv iForveille et al.ll2009H . 

GJ 176 joins a growing number of systems with 
RV-detec ted exoplanets and long-per i od magnetic cy¬ 
cles (e .g. iDumusaue et al.l I2011L 120121 : iRobertson et al.l 
l2013bH . In general, activity cycles which create RV 
signals-presumably via changing magnetic inhibition of 
convection-app ear to be common across a wide range of 
spectral types (iLovis et ah! 120111 : iGomes da Silva et al.l 
l2012t IRobertson et al.ll2013al) . These cycles may either 
mimic exoplanets, or significantly alter the measured or¬ 
bital properties of real planets, so it is important to prop¬ 
erly diagnose and correct them when characterizing RV 
systems. 

The rotation-induced RV signals of GJ 176 again il¬ 
lustrate both the difficulty and importance of care¬ 
ful activity analysis and correction when attempt¬ 
ing to identify low-mass planets in the habitable 
zones of M dw arfs. Recent an a lyses of the quiet M 
stars GJ 581 (IRobertson et al.l l2(314ll and GJ 667C 
(IRobertson fc Ma hadevaril l2014f) reinforced the predic- 
tion of iBoisse et al.l (|2011ll that activity-induced RV sig¬ 
nals appear at the stellar rotation period and its har¬ 
monics. For slowly rotating {Prot lOOd) M stars, these 
rotation harmonics are coincident with the periods of 
planets in the habitable zone. For GJ 176, where the rel¬ 
atively faster rotation period shifts most of its harmonics 
inward of the HZ, we see additional activity signals at in¬ 
termediate periods. If these signals also appear in RV at 
the 0.5-1 ms“^ amplitudes expected for both activity sig¬ 
nals and super-Earths in the HZ, it will again create con¬ 
fusion when searching for planets of astrobiological inter¬ 
est. Indeed, the 39-day rotation period and the 112-day 
activity signal roughly encl ose the optimistic habit able 
zone of GJ 176 according to iKopparapu et ^ (|2013ll . It 
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Fig. 9. — a. RV as versus /d for our HET/HRS RVs of GJ 176. The red line gives our linear least squares fit to the relation, b. 
Generalized Lomb-Scargle periodogram of the HRS+HIRES RVs of GJ 176 before (red) and after (blue) subtracting Equation[T] Note: for 
each of these plots, the secular acceleration has been subtracted from the velocities. 


is therefore clearly possible for the period space of an 
M star’s habitable zone to be completely contaminated 
with stellar activity signals, even when the star rotates 
quickly enough that its harmonics are at shorter periods. 

7. SUMMARY 

We have presented a long-term study of magnetic ac¬ 
tivity of the M dwarf GJ 176 based on spectral activity 
indicators and optical photometry. In addition to the 
stellar rotation period and its first harmonic, we identify 
two intermediate-period activity signals at periods close 
to 2 and 3 times the stellar rotation, and a linear trend 
indicative of a long-period magnetic cycle. The magnetic 
cycle and the stellar rotation appear in RV, although the 
rotation signal is only detectable in the HARPS RVs, 
which cover only a small fraction of the time baseline 
provided by our data and the HIRES RVs. Our results 
lend additional confirmation that the sodium resonance 
lines are a powerful tool for identifying and character¬ 
izing activity-induced RV shifts in M stars, as they are 
sensitive to the magnetic cycle and the 70-day periodic¬ 
ity, and correlate with RV, whereas Ha does not. 
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